AMENDMENT 

In The Claims 

Please amend the claims as follows: 



1. (Currently amended) An apparatus for computing multiple integral of a 
multidimensional integrand function A to be integrated with using a vector map f with 
unbounded support which converts an m (m>l)-dimensional vector having real number 
components into an another m-dimensional vector having real number components, by which a 
multidimensional density function p for the a limiting distribution resulting from repeatedly 
applying the map f to an m-dimensional vector u is analytically solvable, said apparatus 
comprising: 

a first storage unit which stores an m-dimensional vector x = (xi, X2, . . . ,x m ); 

a second storage unit which stores a scalar value w; 
\ a first computing unit which computes a vector x' = f(x) = (x'i, x'2, . . . , x' m ,) resulting 

from applying said vector map f to said vector x being stored in said first storage unit; 

a second computing unit which computes a scalar value w 1 = A(x)/p(x) based on said 
vector x being stored in said first storage unit and said scalar value w being stored in said second 
storage unit; 

an update unit which updates the value stored in said first storage unit by storing said 
vector x f computed by said first computing unit on said first storage unit, and updates the value 
in said second storage unit by adding said scalar value w' computed by said second computing 
unit to a value to be stored in said second storage unit; and 

an output unit which computes , when the number of update times by said update unit 

G 

becomes c (c>lV a scalar value s = w/(c+l) based on said scalar value w being stored in said 

second storage unit when the number of update times by said update unit becomes c (c>l), and 

outputs said scalar value s as a result of the multiple integral. 
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2. (Original) The apparatus according to claim 1, wherein said scalar value stored in 
said second storage unit first is a result from dividing a value resulting from applying said 
function A to said m-dimensional vector stored in said first storage unit first by a value resulting 
from applying said density function p to said m-dimensional vector stored in said first storage 
unit first. 

3. (Original) The apparatus according to claim 1, wherein said scalar value stored in 
said second storage unit first is 0; and 

said output unit computes a scalar value s 1 = w/c, and outputs said scalar value s' as the 
result of the multiple integral instead of said scalar value s. 

4. (Original) The apparatus according to claim 1 further comprising: 

a convergence rate obtainer which obtains convergence rate of scalar values sequentially 
output by said output unit while varying said number of update times c for each of plural vector 
maps gi, g2, . . . , gk (k>2) which are prepared as said vector map f; 

a vector map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects a vector map gh (l<h<k) which shows fastest convergence rate; and 

an output controller which controls said output unit to output said scalar value with using 
said vector map gh as said vector map f and the number of update times c' (c'>c) instead of said 
number of update times c. 

5. (Currently amended) The apparatus according to claim 1, wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 

u, f(u), f(f(u)), f(f(f(u)))), ... 
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resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (ui, U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p(u)= nr 0 m Pi(uO; 
Pi(Ui) ~ c_i|ui|" (1+a) fcrUj— ► — oo; 
pi(Ui) - c_i|uij" (1+a) forUi-> + oo; 
(a>0 5 l<i<m, c_j>0, c +i >0X 



6. (Currently amended) The apparatus according to claim 5 5 wherein said vector 
map f is defined as 

f(u)=(fi(ui), f 2 (u s ), . . . , f m (u m )) 

by a function fj(t) = gi(di t)/d (dj>0) which is defined in l<i<m, and said map gi is 
defined by any one of the following maps (pj(l<j<8) and a natural number n x (ni>2), as follows: 



gi(<Pj(e))=<Pj(rii(e); 




-sgn(tan0)/|tan0| l/a ; 


92(9) = 


-sgn(tan0) x |tan0| 1/a ; 


(p 3 (0) = 


-sgn(cos0) / |tanG| 1/a ; 


94(0) = 


-sgn(cos0) x |tan0| 1/a ; 


9s(0) = 


sgn(cos9) / |tan0| 1/a ; 


96(0) = 


sgn(cos9)x|tan0| 1/a ; 


97(0) = 


sgn(sin0)/|tan9| 1/a ; 


9s(0) = 


sgn(sin9)x|tan0| 1/a ; 


sgn(t) = 


= 1 for t >0; 


sgn(t) = 


= 0 for t =0; 


sgn(t) = 


= -1 for t <0, 



4 




7. (Original) The apparatus according to claim 5 further comprising: 

a convergence rate obtainer which defines said map f for each of plural positive numbers 
qi, q2 9 . . . , q3 (k>2) prepared as an invariable a, and obtains convergence rates of the scalar 
values sequentially output by said output unit while varying said number of update times c; 

a positive number selector which refers the convergence rates obtained by said 
convergence rate obtainer, and selects an integer qh (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said map f with using said positive number qh said 
invariable a, and controls said output unit to output said scalar values with using the number of 
update times c' (c f >c) instead of said number of update times c. 

8. (Original) The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map gi with using plural ones of said 
maps (ft, and obtains convergence rates of the scalar values sequentially output by said output 
unit while varying said number of update times c; 

a map selector which refers to the convergence rates obtained by said convergence rate 
obtainer, and selects one of said maps q)j which shows the fastest convergence rate; and 

an output controller which defines said map gi with using said map cpj selected by said 
map selector, and controls said output unit to output the scalar values with using the number of 
update times c f (cVc) instead of said number of update times c. 

9. (Original) The apparatus according to claim 6 further comprising: 

a convergence rate obtainer which defines said map gi relating to each of plural natural 
numbers Pi, P2, . . . , Pk (k>2) as said natural numbers n^ and obtains convergence rates of the 
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scalar values sequentially output by said output unit while varying said number of update 
times c; 

a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number Ph (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said natural number map gi with using said natural 
number Ph as said natural number n\, and controls said output unit to output the scalar values 
with using the number of update times c 1 (c f >c) instead of said number of update times c. 

10. (Original) The apparatus according to claim 1, wherein said output unit computes 
said scalar value s each time said update unit carries out update, compares said latest scalar value 
s with the former scalar value which is computed at former update, and outputs said latest scalar 
value s if a result of the comparison satisfies a predetermined condition for terminating the 
computation. 

11. (Currently amended) A method for computing multiple integral of a 
multidimensional integrand function A to be integrated with using a vector map f with 
unbounded support which converts an m (m>l) -dimensional vector having real number 
components into an another m-dimensional vector having real number components, by which a 
multidimensional density function p for the a limiting distribution resulting from repeatedly 
applying the map f to an m-dimensional vector a is analytically solvable, said method using 

a first storage unit which stores an m-dimensional vector x = (xi, X2, . . . , x m ), and 
a second storage unit which stores a scalar value w, 
said method comprising the steps of: 
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computing a vector x' = f(x) = (x'i, x f 2, . . . , x' m ) resulting from applying said vector map 
f to said vector x being stored in said first storage unit; 

computing a scalar value w' = A(x)/p(x) based on said vector x being stored in said first 
storage unit and said scalar value w being stored in said second storage unit; 

updating the value stored in said first storage unit by storing said vector x computed by 
^said-first computing unit on said first storage unit, and updating the value in.^aid second storage 
unit by adding said scalar value w computed by said second computing unit to a value to be 
stored on said second storage unit; and 

computing , when the number of update times by said update unit becomes c (c>D, a 
scalar value s = w/(c+l) based on said scalar value w being stored in said second storage unit 
when the number of update times by said update unit becomes c (c>l), and outputting said scalar 
value s as a result of the multiple integral. 



12. (Original) The method according to claim 1 1, wherein said scalar value stored in 
said second storage unit first is a result from dividing a value resulting from applying said 
function A to said m-dimensional vector stored in said first storage unit first by a value resulting 
from applying said density function p to said m-dimensional vector stored in said first storage 
unit first. 

13. (Original) The method according to claim 1 1, wherein said scalar value stored in 
said second storage unit first is 0 ? and 

said outputting step computes a scalar value s' = w/c 5 and outputs said scalar value s 1 as 
the result of the multiple integral instead of said scalar value s. 
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14. (Original) The method according to claim 1 1 further comprising the steps of: 
obtaining convergence rate of scalar values sequentially output by said outputting step 

while varying said number of update times c for each of plural vector maps gi, g2, . . . , (k>2) 
which are prepared as said vector map f; and 

referring to the convergence rates obtained by said convergence rate obtainer, and 
selecting a vector map gh (l<h<k) which shows fastest convergence rate, and 

said outputting step outputs said scalar value with using said vector map gh as said vector 
map f and the number of update times c* (c'>c) instead of said number of update times c. 

15. (Currently amended) The method according to claim 11, wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 

resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (ui, U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p (u) = n i= o m Pi(uO; 
Pi fa) ~ c_i |uj|" (1+a) for Uj -> - oo); 
Pi (Ui) ~ c_j |ui|" (1+a) for Ui -» + oo); 
(a>0, l<i<m, c-j>0, c +i >0). 

16. (Currently amended) The method according to claim 15, wherein said vector map 
f is defined as 

f(u) = (fl(U|), f 2 (u s ), . . . , f m (u m )) 

by a function fj(t) = gj(dj t)/dj (di>0)) which is defined in l<i<m, and said map gi is 
defined by any one of the following maps <f>j (l<j<8) and a natural number nj, (nj> 2) as follows: 
gi (pj(0)) = (pj(ni9); 
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• 


9i(9) = 


-sgn(tane)/|tan0| 1/a ; 


92(6) = 


-sgn(tanG) x |tan9| 1/a ; 


93(6) = 


-sgn(cose)/|tan0| 1/a ; 


94(9) = 


-sgn(cos6) x |tan9| l/a ; 


9s(6) = 


sgn(cos9)/|tan9| 1/a ; 


96(6) = 


sgn(cos9)x|tan9| 1/a ; 


97(6) = 


sgn(sin9)/|tan9| 1/a ; 


9s(6) = 


1 In. 

sgn(sin9)x|tan9| 1/a ; 


sgn(t) = 


= 1 for t >0; 


sgn(t) = 


= 0 for t =0; 


sgn(t) = 


= -1 for t <0 X 



17. (Original) The method according to claim 15 further comprising the step of: 
defining said map f for each of plural positive numbers qi ? q2 ? . . . , qk (k>2) prepared as 

an invariable a, and obtaining convergence rates of the scalar values sequentially output by said 
output unit while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a positive number qh (l<h<k) 
which shows the fastest convergence rate; and 

defining said map f with using said positive number qh as said invariable a, and 
controlling output of said scalar values with using the number of update times c f (c'>c) instead of 
said number of update times c. 

18. (Original) The method according to claim 16 further comprising the steps of: 
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defining said. map gj with using plural ones of said maps cpj, and obtaining convergence 
rates of the scalar values sequentially output by said output step while varying said number of 
update times c; 

referring to the obtained convergence rates, and selecting one of said maps q)j which 
shows the fastest convergence rate; and 

defining said map gi with using said selected map q>j selected, and controlling output of 
the scalar values with using the number of update times c* (c'>c) instead of said number of update 
times c. 

19. (Original) The method according to claim 16 further comprising the steps of: 
defining said map gi relating to each of plural natural numbers Pi, P2, . . . , Pk (k>2) as 

said natural numbers n h and obtaining convergence rates of the scalar values sequentially output 
by said outputting step while varying said number of update times c; 

referring to the obtained convergence rates, and selecting a natural number Ph (l<h<k) 
which shows the fastest convergence rate; and 

defining said natural number map gi with using said natural number Ph, as said natural 
number nj, and controlling output of the scalar values with using the number of update times 
c' (c'>c) instead of said number of update times c. 

20. (Original) The method according to claim 11, wherein said outputting step 
computes said scalar value s each time said updating step caries out update, compares said latest 
scalar value s with the former scalar value which is computed at former update, and outputs said 
latest scalar value s if a result of the comparison satisfies a predetermined condition for 
terminating the computation. 
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2 1 . (Currently amended) A computer readable recording medium storing a program 
for computing multiple integral of a multidimensional integrand function A to be integrated with 
using a vector map f with unbounded support which converts an m (m>l) -dimensional vector 
having real number components into an another m-dimensional vector having real number 
components, by which a multidimensional density function p for the a limiting distribution 
resulting from repeatedly applying the map f to an m-dimensional vector u is analytically 
solvable, said program comprising instructions for caus e s a comput e r to function as : 

a first storag e unit which stor e s storing an m-dimensional vector x = (xj, X2, . . . , x m ) in a 
first storage unit ; 

a s e cond storag e unit which stor e s storing a scalar value w in a second storage unit ; 

a first computing unit which comput e s computing a vector x' = f(x) = (x'i, x f 2, . . . x' m X 
using a first computing unit, the vector x' resulting from applying said vector map f to said vector 
x being stored in said first storage unit; 

a s e cond computing unit which comput e s computing a scalar value w ! = A(x)/p(x) , using 
a first computing unit, the scalar value w' based on said vector x being stored in said first storage 
unit and said scalar value w being stored in said second storage unit; 

an updat e unit which updat e s updating the value stored in said first storage unit , using an 
update unit, by storing said vector x' computed by said first computing unit on said first storage 
unit, and updates the value in said second storage unit by adding said scalar value w' computed 
by said second computing unit to a value to be stored on said second storage unit; and 

an output unit which comput e s computing a scalar value s = w/(c+l \ when the number of 
update times by said update unit becomes c (c>D, the scalar value s based on said scalar value w 
being stored in said second storage unit when the number of update times by said update unit 
becomes c (c>l), and outputs said scalar values as a result of the multiple integral. 
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22. (Original) The recording medium according to claim 21, wherein said scalar value 
stored in said second storage unit first is a result from dividing a value resulting from applying 
said function A to said m-dimensional vector stored in said first storage unit first by a value 
resulting from applying said density function p to said m-dimensional vector stored in said first 
storage unit first. 

23. (Original) The recording medium according to claim 21, wherein said scalar value 
stored in said second storage unit first is 0; and 

said output unit computes a scalar value s' = w/c 5 and outputs said scalar value s* as the 
result of the multiple integral instead of said scalar value s. 

24. (Original) The recording medium according to claim 21, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which obtains convergence rate of scalar values sequentially 
output by said output unit while varying said number of update times e for each of plural vector 
maps gi, g 2 , . . . , gk (k>2) which are prepared as said vector map t 

a vector map selector which refers to the convergence rates obtained by said convergence 
rate obtainer, and selects a vector map gh (l<h<k) which shows fastest convergence rate; and 

an output controller which controls said output unit to output said scalar value with using 
said vector map gh, as said vector map f and the number of update times c* (c'>c) instead of said 
number of update times c. 

25. (Currently amended) The recording medium according to claim 21, wherein a 
multidimensional density function p representing the limiting distribution of a vector sequence 

u,f(u),f(f(u)),f(f(f(u))), 
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resulting from applying said vector map f to a predetermined m-dimensional vector 
u = (uu U2, . . . , u m ) for equal to or greater than 0 times, satisfies the following property of: 
p (u) = n i= o m Pi(ui); 
Pi (uD ~ c_j |ui|" (1+a) for Ui -* - oo); 
Pi (us) ~ c_i |uif (1+a) for u» -> + oo); 
(a>0, l<i<m, c-i>0, c+j>0). 



26. (Currently amended) The recording medium according to claim 25, wherein said 
vector map f is defined as: 

f(u) = (fl(Ui), f 2 (u s ), . . . , fm(Um)) 

by a function f|(t) = gj(dj t)/di (di>0) which is defined in l^Km, and said map gi is 
defined by any one of the following maps q)j (l<j<8) and a natural number n\ (nj>2) as follows: 



gi <pj(9)) = 9j(niO); 


9.(0) = 


-sgn(tanO) / |tan0| 1/a ; 


q>2(6) = 


-sgn(tan0) x |tan0| 1/a ; 


93(6) = 


-sgn(cosO)/|tan0| l/a ; 


94(0) = 


-sgn(cos0)x|tanO| 1/a ; 


9s(0) - 


sgn(cos0)/|tanO| ,/a ; 


9e(0) = 


sgn(cos0) x |tan0| ,/a ; 


97(0) = 


sgn(sinO) / |tanO| 1/a ; 


9s(0) = 


sgn(sinO)x|tan0| 1/a ; 


sgn(t) = 


= 1 for t >0; 


sgn(t) = 


= 0 for t =0; 


sgn(t) = 


= -1 for t <0. 
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27. (Original) The recording medium according to claim 25, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which defines said map f for each of plural positive numbers 
qi, q2, . . . , qk (k>2) prepared as an invariable a, and obtains convergence rates of the scalar 
values sequentially output by said output unit while varying said number of update times c; 

an integer selector which refers the convergence rates obtained by said convergence rate 
obtainer, and selects a positive number qh (l<h<k) which shows the fastest convergence rate; and 

an output controller which defines said map f with using said positive number qh, as said 
invariable a, and controls said output unit to output said scalar values with using the number of 
update times c' (c'>c) instead of said number of update times c. 

28. (Original) The recording medium according to claim 26, wherein said program 
further causes said computer to function as: 

a convergence rate obtainer which defines said map gj with using plural ones of said 
maps q)j and obtains convergence rates of the scalar values sequentially output by said output unit 
while varying said number of update times c; 

a map selector which refers to the convergence rates obtained by said convergence rate 
obtainer, and selects one of said maps q)j which shows the fastest convergence rate; and 

an output controller which defines said map gi, with using said map (pj selected by said 
map selector, and controls said output unit to output the scalar values with using the number of 
update times c' (c'>c) instead of said number of update times c. 

29. (Original) The recording medium according to claim 26, wherein said program 
further causes said computer to function as: 
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a convergence rate obtainer which defines said map gl relating to each of plural natural 
numbers Pi, P2, . . . , Pk (k>2) as said natural numbers nj, and obtains convergence rates of the 
scalar values sequentially output by said output unit while varying said number of update 
times c; 

a natural number selector which refers to the convergence rates obtained by said 
convergence rate obtainer, and selects a natural number Ph (l<h<k) which shows the fastest 
convergence rate; and 

an output controller which defines said natural number map gi with using said natural 
number P h as said natural number ni, and controls said output unit to output the scalar values 
with using the number of update times c* (c'>c) instead of said number of update times c. 

30. (Original) The recording medium according to claim 21, wherein said output unit 
computes said scalar value s each time said update unit caries out update, compares said latest 
scalar value s with the former scalar value which is computed at former update, and outputs said 
latest scalar value s if a result of the comparison satisfies a predetermined condition for 
terminating the computation. 



3 1 . (Original) The recording medium according to claim 2 1 , wherein said recording 
medium is a compact disc, a floppy disk, a hard disk, a magneto-optical disk, a digital versatile 
disc, a magnetic tape, or a semiconductor memory. 
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